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We investigate the paramagnetic periodic Anderson model using the dynamical mean-field theory 
in combination with the modified perturbation theory which interpolates between the weak and 
strong coupling limits. For the symmetric PAM, the ground state is always a singlet state. However, 
as function of the hybridization strength, a crossover from collective to local Kondo screening is 
found. Reducing the number of conduction electrons, the local Kondo singlets remain stable. The 
unpaired /-electrons dominate the physics of the system. For very low conduction electron densities, 
a large increase of the effective mass of the quasiparticles is visible, which is interpreted as the 
approach of the Mott-Hubbard transition. 
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I. INTRODUCTION 

The periodic Anderson model (PAM) is one of the stan- 
dard models for heavy fermion systemslj. In its simplest 
form, it describes a system of localized electronic states 
which hybridize with an uncorrelated conduction band. 
Except fc* a few statements concerning ground state 
propertiescl, no exact solution of the PAM is known up 
to now. TkejKJwujpi been many approxim at e calc ulations 
like QMCMfllra, perturbaW 
as other analytical approachesEj , EJIl3 
However, it is still far from being understood. Closely re- 
lated to-|the PAM is the single-impurity Anderson model 
(SIAMjI The SIAM is one of the best-understood mod- 
els in theoretical physics; besides a broad range of ap- 
proximate solutions even some exact calculations are pos- 
sible (for a detailed review, see reference [j]). The main 
result of these calculations is the emergence of a new tem- 
perature scale Tk {Kondo Temperature), which governs 
the low-temperature physics. For T < Tk, the magnetic 
moment of the impurity is screened by conduction elec- 
trons {Kondo screening). All thermodynamic properties 
can be expressed in terms of Tk- 

However, due to the periodicity of the localized states 
in the case of the PAM, new and more complicated physi- 
cal properties will emerge. The most prominent example 
is the so=called exhaustion problem first mentioned by 
Nozierestil His argumentation was based on the Kondo 
lattice model (KLM), which can be derived from the 
PAM under the condition of half-filling for the legalized 
states using the Schrieffer- Wolff transformationliO. In 
the KLM, where the charge degrees of freedom of the lo- 
calized states have been removed, Kondo screening man- 
ifests itself in the quenching of the magnetic moment of 
the localized spins by the conduction electrons. Reduc- 
ing the number of conduction electrons below half-filling, 
however, the situation changes. The conduction band 
occupation is not sufficient to screen the localized spins 
completely. The natuxe of the ground state is not clear. 
As noted by Noziereso, the situation is different for the 
case of small and large J, respectively. Whereas for small 
J, the screening is a collective effect, in the limit J — > oo, 
it can be understood as the formation of local Kondo sin- 
glets, in which one conduction electron and one S = \- 



spin form a bound singlet state. On removing conduc- 
tion electrons, some of these local Kondo singlets (LKS) 
will be broken. The remaining "bachelor" spins can be 
described as a system of hard-core fermions, for which 
double occupancy is forbiddentJ. The small- J limit with 
its collective screening is more complicated. But nev- 
ertheless, following Nozieres, this case can, similarly be 
mapped onto an effective Hubbard modeled. 



Contrary to the above-described KLM, the localized 
states also have a charge degree of freedom in the PAM, 
and its physical properties are therefore more compli- 
cated. The Schrieffer- Wolff transformation is valid only 
in the small J = limit (V is the hybridization 

strength between localized and the conduction band 
states and U is the on-site Coulomb interaction among 
the localized electrons (see section II). So, especially 
in the limit of strong hybridization between localized 
and conduction states, the KLM is not a priori justi- 
fied as an effective model for the PAM. For the weakly 
hybridized case, the exhaustion in.ijhe.EAM has recently 
been discussed by several authorsuu'lIj'El They find that 
Nozieres' picture of the effective Hubbard model also 
holds in this case and explains the emergence of a new 
low-temperature scale. 



In this paper, we want to present a systematic in- 
vestigation of the PAM. In the following chapter, we 
will introduce the theoretical approach, which fcrni-, 
bines the dynamical mean-field theory (DMFT)E30'E!l 
with a-jwd c.ation of the iterated perturbation theory 



(IPTJ 



tion [II 



OEj. The results will be presented in see- 
In the first part of the discussion, we will 
focus on the hybridization-strength dependence of the 
symmetric PAM, and in the second part, we reduce the 
number of conduction electrons and investigate the ex- 
haustion problem, thereby extending the discussion of 
reference |l9| 
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II. THEORY 



The PAM is defined by its Hamiltonian 

H = E <%*L s s,° + E e /4./^ + a) 
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with s^ )CT (/i.cr) and sj CT being the conduction 



band (/-level) electron annihilation and creation oper- 
ators {nYl = fj^fi^)- e (fc) is the dispersion of a non- 
degenerate s-type conduction band and e/ denotes the 
position of the /-level with respect to the center of grav- 
ity of the conduction band. The hybridization V is taken 
as a fc-independent constant, and finally U is the local 
Coulomb interaction between two /-electrons at the same 
lattice site. 

To obtain the the /-electron Green function G\H (E) — 
a two-step procedure. The 
as dynamical mean-field theory 



(ifialft)), we apply 
first of these— known 
(DMFT)cjO'E3, is a mapping of the PAM onto a simpler 
model, namely the single-impurity Anderson model. The 
second step of our procedure is to find an approximate 
solution of thp_SIAM using the modified perturbation 
theory (MPT)EEffl. 

The starting point of the DMFT is the assumption of 
a k- independent, i.e. local self-energy E CT (i?). It can be 
shown that in this case, the self-energy of the PAM is 
equivalent to the self-energy of a properly defined im- 
purity model (SIAM). The conduction band within this 
SIAM has to be determined by the following expression 
for the hybridization function!!! 

A a (E) =E-e f - £ CT (£) - (G { £(E)) ' (2) 

In the original SIAM, this function is given as A(E) = 
All information of the conduction band and 



E-e(k) ' 

its hybridization with the impurity is contained in A(_E). 
Therefore its knowledge is sufficient to define the electron 
bath of the SIAM also within the DMFT. 

From the self-energy T, a (E) of the SIAM defined by 
A a (E), the PAM /-electron Green function can directly 
be obtained as 



G ul( E ) = jyE 



E-e f - 



v 2 



E-e(k) 



K - V„{E) 



(3) 



Since G\H (E) enters expression (||) , a self-consistent so- 
lution has to be found by iteration. As already noted 
above, the DMFT-procedure becomes exact for a local, 
i. e. ^-independent self-energy. It has been shown that 
for the limit of infinite dimensions, or equivalently in the 
limit of the lattice coordination number going to infin- 
ity, this is indeed exactly the caseEjCJ. Furthermore, 



solving the PAM in three dimensions using perturbation 
theory, it was shown that the results obtained in the local 
approximation compare qualitatively and quantitatively 
very well with-those where the full fc-dependence has been 
consideredtilL-l 

Now the actual problem is shifted to solve the SIAM 
which is defined by A a (E) (see equation (||)). For 
this ta.sk. we. apply the modified perturbation theory 
(MPT)EiEilH which is bas^d, on the following ansatz for 
the electronic self-energy 29 



where (E) denotes the second-order contribu- 

tion to the second-order perturbation, theory around the 
Hartree-Fock solution (SOPT-HFjBEaM Please note 
that the ansatz (^) is ft-independent by construction, the 
basic assumption of the DMFT procedure is therefore al- 
ready incorporated. Hence, the ansatz (Q) together with 
the proposal to fit the parameters as will be described 
in detail below is the only approximation necessary to 
obtain the /-electron Green function. 

The coefficients a a and b a are determined so that the 
first four (?i 6 {0, 1, 2, 3}) moment sum rules 

M(») = fdEE n A^(E) = <[[...[/„, HU. ..,H]_,ft} + ) 

J " v ' 



U(n U l) + 



v (SOC) 



(E) 



b a T,f OC) (E) 



(4) 



-fold commutator 



(5) 



AW(E) = -±!SG£(E + iO+) 



are fulfilled!!!. Qx denotes the imaginary part of x. 

(n) 

Since the moments determine the high-energy 

expansion of the Green function, the compliance of the 
ri = 3 sum rule automatically leads to the correct be- 
havior of G\{l(E) up to the order -grEl. Furthermore, 
the n — 3 sum rule is directly related to the correct po- 
sitions and spectral weights of thts-eharge excitations in 
the strong-coupling limit U — > od23 : E3. This is ensured 
by the occurrence of a higher-order correlation function 
called bandshift B a , which is discussed in detail in the 
context of the Hubbard model in reference |3^, and in 
reference ^0| for the effective SIAM within the dynamical 
mean-field theory. 

Since we use the perturbation theory around the 
Hartree-Fock solution to determine the T, a (E), an- 
other parameter enters the calculation, namely the chem- 
ical potential within the Hartree-Fock calculation. It is 
a priori not evident that this chemical potential should 
be identical to the chemical potential of the full problem. 
In fact, several other choices to determine this parameter 
seem possible. Within the iterated perturbation theory 
(IPT)BEl, which is away from haU-ifilling also based on 
the ansatz (Q), tht; Luttinger sumcIJ or equivalently the 
a rulecic 2 ! is u: 



Friedel sum 



used. This, however, limits the 
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calculation to T = from the very beginning. We there- 
fore define a different constraint, demanding that the im- 
purity occupation number within the Hartree-Fock cal- 
culation is equivalent to the true occupation number. A 
detailed investigation of this choice and other possibili- 
ties is found in reference |3^. Investigating the well-known 
single-impurity Anderson model to test the quality of our 
methodLj, we found that the MPT fulfills the Friedel sum 
rule within numerical accuracy not only under symmetric 
parameter conditions but also in a broad range of param- 
eters, especially when reducing the conduction electron 
density. But contrary to the IPT, the MPT is applicable 
also at finite temperatures. 

Summarizing the features of the MPT, it should be 
considered trustworthy for small U since it is based on 
perturbation theory. But furthermore it is well justified 
in the strong coupling regime, where the main features - 
the charge excitations are correctly reproduced. This is 
clearly one step beyond similar methods which determine 
the parameters, a- and b a with respect to the "atomic" 

limit v = dBBB. 

To calculate the susceptibility, we apply an exter- 
nal magnetic field B cxt which couples to the /- and 
the s-electrons. The susceptibility ^( tot ) is given as 
Qg 1 |_B cxt =o> where M is the total magnetization of the 
system. Since the /- and s-magnetization can be com- 
puted separately, the respective contribution of the f-(s-) 
electrons to ^( tot ) can be determined. 

With the above-described theory, the /-electron Green 
function and all quantities deriveable from it can be 
obtained. This includes several two-particle correla- 
tion functions as e.g. the /-electron double occupancy 
(ni^n^2)0@. However, for the discussion below, we 
will also be interested in two-particle correlation func- 
tions which are not readily obtained from G ii( ^(E). To 
determine these, as e. g. the s-f density correlation func- 
tion (nS n^ S / ) we need a further approximation: We con- 
struct the following effective medium Hamiltonian: 

h rr h& 

(6) 

i.a 

with E CT (E) being the fully self-consistent solution of the 
DMFT-MPT scheme. Since the effective medium Hamil- 
tonian is bilinear in fermion operators, all Green func- 
tions of interest can be evaluated exactly. By construc- 
tion, the single-particle properties of model (|6|) are equiv- 
alent to the original model (]l|) solved within the DMFT- 
MPT scheme. Although we are aware of the fact, that 
using the effective Hamiltonian (|^) must be seen as an 
approximation to the original model ([!]) , it is in our opin- 
ion clearly one step beyond standard first-order pertur- 
bation theory for the latter, which would not reproduce 
the non-trivial results we will discuss below. Let us al- 
ready point out here that these non-trivial results will 
always be accompanied by special features in quantities 



which were derived from the full Hamiltonian ([[]) using 
the DMFT-MPT (e. g. the susceptibility or (n^n^l) ). 
They seem therefore to be not effects due to the replace- 
ment of Hamiltonian with (||) . 



III. RESULTS AND DISCUSSION 

A. The symmetric PAM 

The symmetric PAM is defined by complete particle- 
hole symmetry, i. e. e/ = — y and a particle- hole sym- 
metric density of states of the conduction band with the 
chemical potential located at its center of gravity. 
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FIG. 1: s- and /-densities of states for U = 2, e/ = — y = — 1, 
V = 0.2, rv-f' — n' s ' = 1 and various temperatures (solid 
line: /-DOS, dotted line: s-DOS), the temperature scale is 
explained in the text. 



In figure 1, the s- and /-densities of states (s- and /- 
DOS) with C/ = 2 and V = 0.2 are plotted for various 
temperatures. The energy scale is defined by the free, i.e. 
unhybridized conduction band of unit width and semiel- 
liptic shape centered at E = 0. Within this energy scale, 
the temperatures will be given in We have plotted 
the projections onto the /(s)-states using solid (dotted) 
lines. 

The DOS consist of the charge excitations approxi- 
mately located at e/ and e/ + U which are dominantly 
of /-character and the conduction band mostly of s- 
character which is slightly deformed due to the hybridiza- 
tion. For low temperatures, an additional feature appears 
in the vicinity of E = /i = 0, the Kondo resonancetl It is 
split by the coherence gap which originates from the co- 
herent hybridization between /- and s-states at all lattice 
sites. 

This coherence gap might be the theoretical equivalent 
of the experimentally seen "pseudo-gap" e.g. in SrnBJtl 
or in therSp-called Kondo insulators, e.g. Ce^Bi^Pt^A or 
CeNiSetB. It can be understood as a level-repulsion be- 
tween the conduction band states and the effective /-level 
located at the chemical potential. Since this effective /- 
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level is clearly correlation-induced, the coherence gap is 
as well. 
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FIG. 2: The on-site double-occupancy correlation functions 
and the local spin-spin correlation function (S\ §\ ) as func- 
tion of the hybridization strength V for T = (other param- 
eters as in fig. 1) 



The DOS obtained by our method comgare very, well 
with those calculated by QMC in d = 1B0, d = 2B and 
d = ooHLj. At least for d = 1 and d — 2, there is one 
qualitative difference, however. In the MPT, the Kondo 
resonance is of pure /-character, whereas in the cited 
papers, the conduction band also contributes to the res- 
onance—Whether this is due to the maximum-entropy 
methodu3 necessary to complement the QMC formalism, 
or an artefact of our method, remains an open question. 
However, physically relevant is only the total density of 
states (/ plus conduction band). The choice of project- 
ing the DOS onto the /- and s-states, i.e. onto the basis 
given by the V = solution of the problem is rather ar- 
bitrary. Therefore, the above-discussed differences seem 
to be of minor importance. 

Focusing on the Kondo screening problem, the ques- 
tion arises on how one can observe it. One possibility 
found in the literature is the definition of an effective 
magnetic moment T\{T). This is motivated by the Curie 
law which, however, only holds for high temperatures. 
In the SIAM, there is a suppression of T\{T) coinciding 
with the temperature scale Tk which also governs the 
other low-energy properties of the system. In analogy to 
the SIAM, the behaviour of T\(T) is often interpreted as 
indirect indicator for the onset of Kondo screening in the 
PAM0Q. Another criterion of Kondo screening, which 
we will focus on, might be found in spin-spin correlation 
functions. Let us look for example at the problem of two 
S = h spins. Of the four possible states of this system, 
three are of triplet and one is of singlet nature. The spin- 
spin correlation function takes the value (S a Sb) — 4 for 

the triplet states and (S a Sb) = — | for the singlet state 
(a and b denote the two spins). In the following, we will 
discuss only the on-site interband spin-spin correlation 
function (S^'§^). It is obvious that this function can 
be an indicator only for local singlet formation, as will 
be discussed in more detail below. 



In figure 0, the on-site interband spin-spin correlation 
function (S^'Sj 8 ') as well as several on-site double occu- 
pancy correlation functions are plotted as function of the 
hybridization strength V for T = (other parameters as 
in figure 1). With increasing V, the on-site spin-spin cor- 
relation function approaches the value — |. In the same 
range of V, the interband double occupancy with parallel 
spin, {n^Jrifa) vanishes, whereas the respective correla- 
tion function with opposite spin indices, {nfjnf}^) stays 
almost constant. In analogy to the two-spin problem, 
these are clear indications of a local singlet correlation. 
A similar transition in the spin-spin correlation function 
has been seen in a PAM with next-neighbor hybridiza- 
tion using|XJMC, where it has been interpreted as singlet 
formationS. 

In the upper graphs of figures || and |J, the correlation 
functions of figure ^| are plotted as function of tempera- 
ture for large (V = 0.6) and small (V = 0.2) hybridiza- 
tion strengths. Additionally, the lower graphs show the 
respective susceptibilities. In the large- V case, the above 
described transition in the various correlation functions 
is clearly visible: (5^ ) — ► — | and {n\^n[ s J) shows 
a huge drop around T w 3000, being of very small value 
at low temperatures. In the same temperature range, the 
susceptibility vanishes. Both the /- and s-contribution 
to ^( tot ) disappear simultaneously. From these findings, 
we propose the occurrence of local Kondo singlet forma- 
tion. With the term local Kondo singlet (LKS) , we want 
to stress, that the singlet formation is predominantly a 
local process determined by the binding of one /- and one 
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FIG. 3: Upper panel: the correlation functions of figure ^, 
but as function of temperature for V = 0.6. Lower panel: 
the susceptibility ^' to *' as well as the /- and s-contribution 
to X (tot) are plotted. The thin dotted line (x (frcc band) ) shows 
the susceptibility of a free, i.e. unhybridized and uncorrelated 
conduction band with same parameters as the s-band used in 
the other calculations. 
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s-electron at each lattice site. In the opposite case of col- 
lective Kondo screening as it could occur for small V, the 
local correlation functions discussed here need not show 
any particularities. Our proposal is further supported by 
the behavior of the conduction band double occupancy 
{n\ s Jnf^ a ). In the large-y region, where we propose the 
LKS formation, this correlation function is reduced com- 
pared to the small- V case, as can be clearly seen in figure 
|2|. This indicates a tendency towards localization of the 
conduction electrons, which is exactly what one would ex- 
pect in the case of LKS formation. The unique tempera- 
ture scale which we identify with the Kondo temperature 
of Tk ~ lOOOif seems to be very large. This cannot be 
due to the "Hartree-Fock-character" implied by the sim- 
ple effective medium approach of equation @ used to de- 
termine the higher correlation functions, since the same 
temperature scale appears in the susceptibility, which is 
determined from the full Hamiltonian (equation (jlj)) in 
the DMFT-MPT scheme. Whether T K is over-estimated 
due to the mean-field character of the DMFT or the use 
of the MPT, or whether it displays the true behavior of 
the system, would be speculation. However, comparing 
our method with other means of calculation (e.^ .numer- 
ical renormalization group calculations (NRG)cj'E3), the 
DMFT-MPT seems to overestimate energy scales. We 
believe this is connected to the fact that the MPT, as any 
perturbative method, is unable to reproduce the expp=. 
nential temperature scale typical for the Kondo physicso 
A detailed comparison with complementary methods, e. 
g. NRG, is needed to shed further light on this. 



a spin gap. However, for small hybridization the vari- 
ous correlation functions discussed above show none of 
the occurrences which led us to the conclusion of the lo- 
cal Kondo singlet formation in the large-I^ case. The 
susceptibility shows two features. The one at the high 
temperature Thigh ~ 4000^ corresponds to the der- 
ealization of the /-electrons due to thermal excitations. 
As can be seen at the decomposition of X : only the 
/-contribution is responsible for this feature. The con- 
duction electron contribution x^ resembles at and be- 
low Thigh still the susceptibility of a system of uncor- 
related electrons (xfreeband)- The derealization of the 
/-electrons is visible in the increase of {n ia rifj^) for 
T > Thigh- The Kondo screening sets in at the much 
lower Kondo temperature Tk ~ 100K. This corresponds 
to the temperature where the Kondo resonance appears 
in the DOS (see figure |l|). As already mentioned, the lo- 
cal correlation functions show here only a weak signature 
of Kondo screening. The Kondo screening for small V is 
a collective effect, the quantities accessible by our theory 
do not allow more detailed investigation of this state. 

To summarize the discussion of the symmetric PAM, 
we find two different kinds of Kondo singlet formation de- 
pending on the hybridization strength. Where for small 
V the singlet formation involves non-local or collective 
screening of the /-electrons by the conduction band elec- 
trons, in the large V domain, the singlet formation is 
dominantly a local process. At every lattice site, one /- 
and one s-electron couple to form a local Kondo singlet 
(LKS). 
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FIG. 4: As figure |[ but with small V = 0.2 



For small V (figure Q), the situation is completely dif- 
ferent. The systems with small and large V behave sim- 
ilar only for very low and very high temperatures. For 
very high temperatures, a Curie-like behavior is found 
as expected. For T = 0, the susceptibility vanishes also 
in the small- F case. This is compatible with the ex- 
act result by Tsunetsugu et al.0, who proved the singlet 
nature of the ground state of the symmetric PAM in- 
dependent of the size of V and showed the existence of 



B. Exhaustion problem 

One important question concerning the screening be- 
havior is that of exhaustion: What happens if the number 
of conduction electrons is reduced (the number of /- 
electrons is fixed n^' = 1)? This .question was recently 
brought into attention by NozieresO In this section, we 
want to present our results concerning "exhaustion" . 

In figure ^ the /-electron density of states (/-DOS) is 
plotted for T = and various . The main change in 
the /-DOS is the shift of the charge excitations as well as 
the Kondo resonance towards lower energies. The shift of 
the charge excitations is due to an adjustment of e/ which 
is necessary to keep the constraint = 1 with de- 
creasing n^ s \ The Kondo resonance also moves towards 
lower energies as its position is pinned to the chemical 
potential fi. The coherence gap stays in the vicinity of 
fx, too. However, the relative position of the gap, the 
Kondo peak and /i changes. Whereas for — 1, the 
Kondo resonance and the gap are centered symmetrically 
around /x, for rS s ^ < 1 the situation becomes asymmet- 
ric. The chemical potential moves into the lower half 
of the resonance. As a direct consequence, the shape of 
the resonance becomes asymmetric. Thia-atrongly resem- 
bles the behavior found for the SIAM£E3. The relative 
shift of fi and the Kondo resonance has also an impor- 
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FIG. 5: /-densities of states for T = in the vicinity of the 
chemical potential (arrows) . The system parameters are as in 
figure |l| but the conduction band occupation is varied from 
n' s ' = 1 to rv 3 ' = 0.2 as indicated. The inset shows the 
/-DOS over the full energy range 



tant consequence for the coherence gap. For < 1 the 
system is metallic since \i is no longer located within the 
gap. With decreasing n^ s \ the gap moves further away 
from the chemical potential. The imaginary part of the 
self-energy QT, a (E) at the gap increases thereby since it 
shows a "Fermi-liquid" E 2 dependence around E = \i. 
For slightly below unity, the gap is still present in 
form of a "pseudo-gap" . However, for approaching 
zero, the gap closes completely. 

This resembles the behavior found in CezEi^Ptz upon 
doping with lanthanumeZI. The undoped system is an 
insulator with a very small gap believed to be a proto- 
typical Kondo insulator. On doping La, the system be- 
comes a metallic heavy-fermion system. In reference 
this is interpreted as due to the de-construction of lattice 
coherence by disorder. From the behavior of the PAM, 
one might also conclude that simply the change of the 
electron density due to doping would also be sufficient to 
explain this metal-insulator transition. 

In figure |^ a broadening of the upper charge excitation 
is observed for — > 0.2. Here, the upper charge exci- 
tation overlaps with the conduction band which results 
in a stronger hybridization. This could be prevented by 
using a larger value for U. However, the impact of this 
on the results discussed below is negligible. 

Let us now turn to the actual problem of exhaustion. 
In the case of the symmetric PAM (n^ — n^' = 1), 
we have found a singlet ground state with aiLnite spin 
gap which is in agreement with QMC resulteEll and ex- 
act statements concerning the ground stataa. It is not 
clear how the system will react when the number of con- 
duction electrons is reduced while n^' — 1 is kept con- 
stant. One possibility would be a partial or complete 
breakdown of the screening since the conduction-band 
filling is not sufficient to screen the magnetic moments of 
all /-electrons. However, conduction-band mediated /- 
/ correlations could again lead to a singlet ground state 
via intersite singlet correlations. This scenario has been 
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FIG. 6: s-DOS and imaginary part of the hybridization func- 
tion — QAc(E) close to the chemical potential \x for the same 
parameters as in figure |H[ The dotted line indicates the po- 
sition of /i for the respective value of n^ 3 \ The curves for 
7^ 1 have been shifted vertically for clarity. 



confirmed for the two-impurity Kondo problemEa. 

Recently, the exhaustion problem Jjia.s Jpegu, .discussed 
as the origin of a newiaw-energy scaleQ0E3E3o'Ej. In one 
of these publications!^, the iterated perturbation theory 
has been applied, which is very similar to the approxi- 
mations used in this paper. There, the authors found a 
gap close to the chemical potential in the effective hy- 
bridization defined by — < ^A a (E) (see equation (^)). In 
figure ||, — < ^sA (T (E) as well as the conduction band den- 
sity of states (s-DOS) are plotted in the vicinity of /j, for 
various conduction band fillings. The parameters corre- 
spond to the DOS plotted in figure ||. In the s-DOS, the 
coherence gap but no Kondo resonance is visible. The 
latter appears only in the projection onto /-states (/- 
DOS) and is therefore of pure /-character. Again, the 
coherence gap shifts together with \x towards lower en- 
ergies on reducing . Let us focus on — ^sA a (E), the 
upper picture in figure ||. Independently of there ex- 
ists a gap/dip close to the chemical potential. However, 
in the symmetric case (n^ = 1) a sharp 5-like peak ap- 
pears exactly at E = /x. So — ^A a (E = /i) is large in 
the symmetric case, but already for any small change in 
n( s \ —^sA a (E = n) becomes small since fi lies within 
the gap or dip. In reference H9l the authors interpret 
the dip which they find for nP> = 0.4 in the effective 
hybridization as sign of the exhaustion of the conduction 
electrons. However, this reduction of — c iA a (E) around 
E = fx, which we also find for — 0.4, continuously 
develops into a gap for = 0.9. In the interpretation 
of the cited work, this would imply that the exhaustion 
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FIG. 7: The T — O-susceptibility and its respective /-level 
and conduction band contributions for small V = 0.2 (up- 
per panel) and large V — 0.6 (lower panel) as function of 
conduction band fillig. In the inset of the upper panel, the 
effective mass is plotted for both values of V. The remaining 
parameters are as in figure ti. 



In figure 0, we have plotted the zero-temperature sus- 
ceptibility for small V = 0.2 and large V — 0.6 as func- 
tion of the conduction band occupation number n 1 --^ . In 
addition to that, the effective mass m* = 1 
plotted in the inset of the upper panel. 

From these figures, one gains some insight in the sta- 
bility of the LKS in case of exhaustion (n^ < n^' = 1). 
In the small- V case, x^ ^ is roughly proportional to the 
value of the DOS at E = /i for less than half-filled con- 
duction band. Both the /- and s-electrons contribute to 
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Contrary to that, for the LKS-dominated system 



(V = 0.6), only the /-electrons react to an external field. 
The conduction band contribution 1S negligible for 
< 0.9. So we concluse that all conduction electrons 
are more or less bound into LKS states and therefore in- 
able to react to an external magnetic field. Another hint 

can be derived from the fact that (S^S^') decreases lin- 
early with (not plotted). This can be interpreted as 
a linear decrease of the number of LKS. I.e. all available 
conduction electrons form a LKS, and n' bach ) = 1 — 
/-electrons remain unpaired. These are the equivalent of 
the "bacheloij-spins" discussed by Nozieres for the Kondo 
lattice modeEJ. Following his argumentation, they can 
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FIG. 8: The susceptibility, its /- and s-contributions and the 
on-site interband spin-spin correlation function (Sj^S^) as 
function of temperature for three different conduction band 
occupations n' s ' £ {0.9, 0.8, 0.4}. The hybridization strength 
is V = 0.6, the data corresponds to figure (?]. Please note that 
we have plotted the negative of (Sj^ ) and rescaled it by 
a factor 10 for the upper two, and a factor 1000 for the lower 
plot. 



fermions. In 



be described as a system of n^ bach ' spin-i 
the KLM these are "hard-core" -fermions since for the 
localized spins, double occupancy is strictly forbidden. 
In our case of the PAM, double occupancy is princi- 
pally possible, only suppressed by the on-site repulsion 
U. Therefore, these bachelor fermions should show simi- 
larities to a non-degenerate Hubbard model with finite U. 
One featu re , of the Hubbard model is the Mott-Hubbard 
transitionEaO. For an exactly half-filled system with suf- 
ficiently large U, the otherwise metallic system becomes 
insulating due to the electron correlations. Approaching 
the half-filled situation from lower carrier concentrations, 
the effective mass strongly increases, and diverges at half- 
filling. Bearing this in mind, the behavior of the effective 
mass shown in the inset of figure [?] can be understood. 
The most prominent feature is the increase of m* for 
n^ s ' — > 0. The latter, however, implies that the above- 
discussed bachelor-fermion model approaches half-filling, 
since the number of unpaired /-electrons goes to unity. 
In terms of this model, the system is in the proximity of 
the Mott-Hubbard transition. Please note, that in the 
PAM with finite hybridization, n/ s ' = is only possible 
for extreme parameter conditions, namely ej — > —oo and 
U — > oo. Our method is not feasible for these parameters. 

Let us now discuss why the susceptibility is enlarged at 
n( s ' = 0.9. Looking at the partial contributions an d 
at — 0.9, we note that both take on large values. 
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However, while \^ is positive, the conduction band con- 
tribution becomes negative (see lowest panel of figure ||). 
We believe this is due to the proximity of a magnetically 
ordered phase. Since we do not allow for antiferromag- 
netic ordering in our method, we obtain a paramagnetic 
system for the shown parameters. For larger values of 
U, however, we also see an onset of ferromagnetic or- 
der. The critical U c is lower for « 0.9 than for other 
values of n^ s \ So the peak in the susceptibility is in 
our opinion due to the proximity of the ferromagnetic 
phase. The opposite signs of x an( i X result from 
the fact that in the ferromagnetic phase, the conduction 
band will be polarized antiparallel to the /-levels, as will 
be shown in a forthcoming publication. However, using 
quantum Monte Carlo naethods, antiferromagnetic order 
was found for — > lu. In principle, the investigation 
of an antiferromagnetic order is also possible within our 
method and is planned for future studies. The results of 
reference [| suggest that for nS^ = 1 and — ► 1 the 
antiferromagnetic phase will be stable against the ferro- 
magnetic phase. 

Note that also in the weakly hybridized case, this in- 
crease in m* is visible. This is similar to the findings in 
reference and can be understood in a similar fashion 
as in the large-y case (a detailed discussion is given in 
reference |J) . 

In figure |^, the temperature dependence of the suscep- 
tibility and the on-site interband spin-spin correlation 
function is plotted for V = 0.6 and £ {0.4,0.8,0.9}. 
As can be seen, the temperature where the LKS forma- 
tion occurs, is nearly independent of n^ s '. The LKS 
formation should not be confused with the potential 
new low-energy scale discussed in recent publicationstlEJ. 
In our results, the only possible indication for this en- 
ergy scale could be the increase of m* . Complementary, 
methods such as numerical renormalization theoryauUH 
should be used to get more insight at this question. 



good accordance with (numerically) exact 
methods has been found. Being fast and numerically sta- 
ble, all parameter regions of the model can conveniently 
be investigated, this includes T — as well as finite tem- 
peratures. 

The density of states generally consists of the charge 
excitations of the localized level, the conduction band 
structure and, for low temperatures, an additional peak, 
the Kondo resonance. At least for symmetric parameters 
the latter is split by the coherence gap. The Kondo reso- 
nance is ascribed to the phenomenon of Kondo screening 
meaning the quenching of the magnetic moment of the 
localized levels by the conduction electrons. 

In the case of the symmetric PAM, where the /-level 
and the conduction band are both half-filled and the two 
charge excitations lie symmetrically around the chemi- 
cal potential, the ground state is always a singlet with a 
spin gap which is in accordance with the assumption of 
complete Kondo screening. Investigating several on-site 
correlation functions and the susceptibility as function of 
the hybridization strength V , we see a crossover between 
two qualitatively different regions. Whereas for small V, 
the Kondo screening is a collective effect, for interme- 
diate and large V, the screening is a dominantly local 
process. At each lattice site a local Kondo singlet exists. 
It is built up by one conduction- and one /-electron spin. 

On reducing the number of conduction electrons, the 
LKS remain stable. However, due to unavailability of 
s-electrons, a finite number of /-electrons are unpaired 
(bachelo£-.fermions) . .-Following a recent reasoning by 
NozieresEil and othersO, the low-temperature physics of 
the system should be describable by an effective model 
which only regards these bachelor fcrmions. Our results 
are compatible with this proposal. In the susceptibility, 
we found indications of the proximity of magnetically or- 
dered phases. These will be subject of a forthcoming 
paper. 



IV. CONCLUSIONS 

In the present work, we have studied the periodic An- 
derson model (PAM) using the modified perturbation 
theory in the context of the dynamical mean-field the- 
ory. This approach is well motivated both for the large 
and small coupling regime. Furtherppre, applying it to 
the single-impurity Anderson modelE3 and the Hubbard 
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